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ABSTRACT 

Long  term  numerical  integrations  were  performed  with  a  two-level 
model  utilizing  the  quasi-geos trophic  equation  set.   Friction  was  in- 
corporated through  the  vertical  derivative  of  the  eddy  stress  and  heat- 
ing was  applied  as  a  linear  function  of  latitude.   The  long-term 
interactions  between  one  or  two  waves  in  the  zonal  direction  and  the 
mean  flow  were  examined. 

The  stability  of  the  initial  mean  flow  was  investigated  with  lin- 
earized equations  and  the  most  unstable  wave  numbers  were  determined. 
The  nonlinear  equations  were  integrated  using  the  most  unstable  wave 
in  the  initial  conditions.   300  day  forecasts  were  made  with  various 
values  of  the  heating,  |3,  wall  separation,  and  wave  number.   A  second 
wave  number  was  introduced  into  the  experiments  in  order  to  determine 
the  most  likely  wave  number  for  the  fully  evolved  mean  flow.   Constant 
amplitude,  baroclinic,  propagating  disturbances  were  obtained  in  every 
case  where  the  wall  separation  was  4,000  km.   The  experiments  with  the 
wall  separation  of  8,000  km  produced  non-steady  waves  with  fluctuating 
baroclinic  and  barotropic  interactions. 


TABLE  OF  CONTENTS 

I.  INTRODUCTION - — -  6 

II.  THE  MODEL  AND  FORECAST  EQUATIONS  8 

III.  PROCEDURE  AND  GOALS 17 

IV.  VERTICAL  MOTION  EQUATIONS —  19 

V.  COMPUTATION  OF  AVERAGE  NORTH-SOUTH  WIND 20 

VI.  ENERGY  AND  ENERGY  TRANSFORMATION  EQUATIONS  — 21 

VII.  INITIAL  CONDITIONS - 24 

VIII.  RESULTS - 26 

A.  HIERARCHY  OF  UNSTABLE  WAVE  NUMBERS  - 26 

B.  INTERACTIONS  BETWEEN  WAVE  AND  MEAN  FLOW  --- 26 

C.  u  AND  uT  VERSUS  HEATING  FOR  THE  STEADY 

WAVE  EXPERIMENTS 39 

D.  A  BAROCLINIC,  STEADY  WAVE 39 

E.  A  NON-STEADY  WAVE  WITH  BAROTROPIC  EFFECTS —  45 

F.  A  TWO  WAVE  EXPERIMENT - 53 

IX.  CONCLUSIONS 61 

LIST  OF  REFERENCES 64 

INITIAL  DISTRIBUTION  LIST - 66 

FORM  DD  1473 72 


TABLE  OF  SYMBOLS  AND  ABBREVIATIONS 


P  Derivative  of  coriolis  parameter  at  mid-latitude 

y  Eddy  viscosity 

C  Specific  heat  at  constant  pressure 

C  Drag  coefficient 

f  Coriolis  parameter  at  mid-latitude 

g  Gravity 

oi  dp/dt 

R  Gas  constant 

p  Density 

T  Average  temperature  from  standard  atmosphere 

Q  (R2T/p2g)  (dT/dz  +  g/Cp) 


C      -fv2* 


o 

z        Height 

H        R/C 

P 

S        Heating  added  per  unit  mass 
W       Distance  between  the  walls 
T       Eddy  stress 


q        Potential  vorticity 


ACKNOWLEDGEMENTS 

The  author  wishes  to  express  his  deep  appreciation  and  thanks  to 
Dr.  R.  Terry  Williams  whose  interest,  counsel,  and  recommendations 
were  truly  exceptional. 

I  would  further  like  to  thank  Dr.  George  J.  Haltiner  for  his 
assistance  and  suggestions  which  contributed  to  this  study.   A  thank 
you  is  also  offered  to  Frank  H.  Taylor,  portions  of  whose  April  1970 
thesis  results  and  model  were  used  in  this  paper. 


I.   INTRODUCTION 

In  a  classic  experiment,  Phillips  (1956)  performed  a  long-period 
numerical  forecast  with  a  two-level  quasi-geostrophic  model  in  a  zonal - 
ly  periodic  channel  on  a  Beta  plane.   Friction  and  heating  were  includ- 
ed in  his  model  with  the  latter  as  a  linear  function  of  latitude. 
Although  Phillips'  experiment  did  not  reach  statistical  equilibrium, 
he  demonstrated  that  the  principal  features  of  the  general  circulation 
of  the  atmosphere  could  be  reproduced  with  a  simple  numerical  model. 
More  refined  studies  involving  the  primitive  equations  have  been  car- 
ried out  by  Smagorinsky  (1963),  Mintz  (1965),  Manabe  and  Smagorinsky 
(1967),  Manabe  (1969),  and  Kasahara  and  Washington  (1967)  and  (1971). 
These  studies  which  more  accurately  model  the  actual  atmosphere  are 
sufficiently  complicated  so  that  it  is  difficult  to  understand  the 
basic  interactions  between" the  mean  flow  and  the  disturbances. 

For  this  reason,  we  wish  to  re-examine  Phillips'  original  simple 
model.   In  this  study,  Phillips'  model  is  further  simplified  to  repre- 
sent the  x-variation  of  the  disturbance  with  one  or  two  Fourier  compon- 
ents. This  restriction  is  partially  justified  by  the  fact  that  one 
predominate  wave  number  was  observed.   Similarily,  in  dish-pan  experi- 
ments with  inner  cores,  single  wave  numbers  were  observed  in  many  in- 
stances; see  Fultz  (1959)  and  Hide  (1953).   Herring  (1963)  and  (1964) 
has  used  the  technique  of  examining  the  interaction  between  a  disturb- 
ance of  sinusoidal  variation  and  a  mean  field  in  his  study  of  convection 
between  two  rigid  plates.   His  estimates  of  the  heat  fluxes  between  the 
plates  were  within  20$  of  the  value  obtained  from  laboratory  experiments. 


This  suggests  that  a  simplified  dynamic  model  of  the  atmosphere  which 
involves  the  interaction  between  a  simplified  disturbance  and  a  mean 
flow  can  give  reasonable  approximations  of  the  northward  fluxes  of  heat 
and  momentum. 

Therefore,  in  this  study,  the  complexities  have  been  kept  to  a 
minimum  to  simplify  interpretation  of  the  interactions.   The  two-level 
model  utilizes  the  quasi-geostrophic  equation  set.   Fluid  motion  in  the 
model  is  restricted  to  the  mean  flow  and  one  or  two  wave  numbers.   The 
mean  flow  and  Fourier  amplitudes  of  the  two  disturbances  are  allowed  to 
vary  in  time  and  y-space;  hence  an  adequate  description  of  barotropic 
as  well  as  baroclinic  interactions  are  possible. 

However,  such  a  drastically  simplified  model  must  omit  many  details 
of  the  real  atmosphere.  In  the  case  of  this  model,  a  number  of  physical 
properties  of  the  earth's  atmosphere  are  missing.  Two  waves  in  juxta- 
position will  not  interact  directly,  but  can  only  influence  each  other 
through  the  modification  of  the  mean  flow.  No  oceans,  moisture,  moun- 
tains, or  east-west  heating  gradient  are  allowed  in  the  model.  Conse- 
quently, this  model  appears  to  be  one  of  the  simplest  models  which  can 
describe  some  of  the  principal  mechanisms  of  the  general  circulation. 


II.   THE  MODEL  AND  FORECAST  EQUATIONS 

The  model  used  in  these  experiments  utilizes  the  simple  two-level, 
quasi-geostrophic  equation  set.   The  atmosphere  is  divided  into  four 
layers  of  thickness  Ap/2,  numbered  0  to  4  in  Figure  1. 
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0 

<»> 

0 

co=o 

4P  — 2 

3 
24p  : = 4 

,    Fig.    1.       Two^level_model .      

The  earth's   surface  is   assumed  to  be  a  smooth,    flat,    land  surface 

with  no  oceans.      Therefore,    the  vertical  motion  at  the  top  and  "bottom 

of  the  atmosphere  is   taken  to  be  zero. 

Internal   friction  is  represented  by  the  vertical  derivative  of  the 
eddy  stress    T  . 

The  quasi-geostrophic  vorticity  and   thermodynamic   equations   are 


!t-+W.VC+Pov-fo^=-lk.Vxg     U 


(2.1) 


dt     dp  dp  p  (2.2) 


KS 

P 


where  H  =  R/C   ,  S  =  net  heating  from  radiation,  and  a,  the  static 
stability,  is 

°*   *P  trd  -  ,1  -  Sp  [§-  +  §I  =  ^VcT*      (2.3) 

P  8  P  g    P 

T  is  taken  from  a  standard  atmosphere. 


Apply  Equation  (2.1)  at  levels  1  and  3,  which  yields 

ST  +  wi  •  vq  +  P0vx  -  £Q  J    "  $,k  'v  *  <V  V        •      <2-4> 


SC3  <n 


^  +  W     •  VC,  +  Pv,  +£    tJ     -  £|k    •  V  x(TA    -    T0)  (2.5) 


If  we  apply  Equation   (2.2)    at   level   2,   we  obtain 


y£>.  +  ^  ■  v^- 


-  KS, 


>p'2    2   vvdpy2    2      Ap  (2.6) 

The  heating  function,  S_,  is  restricted  to  a  linear  function  in  y 

as  utilized  by  Phillips  (1956).   The  model  has  heating  in  the  south  and 
cooling  in  the  north. 

S2  s  h(|  -  y)/w  (2.7) 

where  w  is  the  north-south  wall  separation  distance. 

The  amount  of  net  heating  is  proportional  to  h  which  varies  from 

0.001  to  0.016  kilojoules  ton"  sec"  . 

The  eddy  stress  T  at  the  top  of  the  atmosphere  is  taken  to  be  zero, 
The  eddy  stress  term  T  at  level  2,  takes  the  form 

2     aw 

T2  =  "  P  &   §F  (2.8) 

where  y  is  the  eddy  viscosity.  If  we  follow  Smagorinsky  (1963),  the 
finite  difference  form  becomes 

T2  =  -<P7)2  (gP2/Ap)  (W3  -  Wx)  (2.9) 

where  (p«g/Ap)  =  1/7.9  km   is  the  inverse  depth  of  the  750-250mb  layer 

'  5°      -4      -1  -1 

x  10   ton  m   s   is  an  exchange  coefficient.   The 


and  (pr)2  = 


225 
500 
smallest  estimates  of  (p7)~  was  proposed  by  Rossby  and  Montgomery  (1935) 

for  stable  conditions,  while  the  largest  was  utilized  by  Riehl  (1951). 

The  intermediate  value,  225,  was  proposed  by  Palmen  (1955). 


T2  "  "  C2   <W3  "  V 

(  50 

where  C     = 

J  225 

x  10"    /79   ton  s"1 

(500 

Thus,  the  stress  term  at  level  2  becomes 

(2.10) 

(2.11) 

The  majority  of  the  experiments  were  conducted  using  Palmen's  exchange 
coefficient  of  225, 

The  eddy  stress  term  T,  at  the  surface,  level  4,  as  formulated  by 
Phillips  (1956)  is  assumed  to  be  in  the  direction  of  the  surface  wind 
and  proportional  to  the  square  of  the  wind  speed. 

T4  S|  CD  KK  (2'12) 

where  C  is  the  drag  coefficient  =  0.015. 

The  surface  geostrophic  vorticity  is  approximated  by  the  750mb 
vorticity 

k  -'  £3  -  -r  y2  *3 

o 

thus  W,  =  \V0  . 
4    3 

So,   by  defining  C^  =     |     CD    I \V^  |  (2.13) 

T4  =  C4  \V3  (2.14) 

The  average  value  of  the  750mb  wind  speed  was  taken  to  be  5ms 
Thus,  our  assumptions  make  T,  a  linear  function  of  wind  speed  with 

C4  -  3.69  x  10"5  ton  s"1  (2.15) 

Therefore,  as  boundary  conditions  at  the  top,  internal,  and  bottom 
layers  of  the  atmosphere 

CD   =  GO,  =  0 

o    4 

T   =  0 
o 

T2  -  -  c2  <w3  -  V 
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\  -   C4  W3 


For  convenience,  the  following  quantities  are  defined 
A      *1  +  *3 


♦l  +  *3 
*j  s      2  (2.17) 


implies 


♦l  "   i  +  *T  *3  =  K  '    *T 

Cx  +  C3  =  Ql  -  C3 

tM       2  ^T      2 

Cj  s  bM  +  £T  C3  =  ^  -  CT 

the  mean  quantities  may  be  interpreted  as  applying  at  level  2. 

Integrate  the  hydrostatic  equation  T^  =  "  a  from  level  1  to  3, 
which  yields 

^  =  f-   In  (P3/P1)  (2.18) 

Equation  (2.18)  shows  that  <|L  is  proportional  to  the  mean  temperature 
of  the  layer. 

If  (2.16)  and  (2.17)  are  introduced  into  the  thermodynamic  equation 
(2.6) ,  we  obtain 

ft  *r  +  Vv*r  -  -V  "  J1  <2-19> 

The  mean  vorticity  equation  is  obtained  by  adding  Equation  (2.4) 
and  (2.5)  and  dividing  by  two. 

It  £m+  %  •  VV  \  •  v«t  +  povm  -  §V|k ' v  x  T4  (2-20) 

The  thermal  vorticity  equation  is  obtained  by  subtracting  Equation 
(2.4)  from  (2.5)  and  dividing  by  two. 
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5i«T  +  \  •   ^T  +  WT  '  %I  +  P0VT  -  £„  4  - 

-  §£lk  •  V  x  T2  +  f^lk  ■  V  x  T4        (2.21) 

Solve  (2.19)  for  ao_  and  apply  this  solution  to  the  thermal  vorti- 

city  equation  (2.21). 

Substitute  the  following  geostrophic  relations 

W  =  |-lk  xV(),        C=f"V2^) 
o  o 

into  Equation  (2.20)  and  multiply  by  f  ,  which  becomes  the  Mean  Vorti- 


oTV\  +  f~   (,kxVi>    "  V<V\>   +I   dkxV(|)T)-V(V2«tlr) 


city  Equation: 

9  „2 

o 

+  po§r  "K4  (v2*t  -v2i>  =0  (2-22) 

Introduce  the  geostrophic  relations  for  W  and  Q   into  (2.21)  and 
multiply  by  f  ,  which  gives  the  Thermal  Vorticity  Equation: 

|^(v2-  v2)    (tj,  +  \     (Ik  x  V^)  •  V(V2  -  \12)    ^ 
o 

+  \  (Ik  x  v<|)T)  •  v<y\)   -  Pq  ^  +  K^(v2^  .  v2^} 

o 

2  ^2         2x 
+  M  2"^  =  -  K2  V^  (2.23) 


where 


2f  2 
H2  =  — ~  =  3.48  x  10"12  m"2  (2.24) 

(Ap)  C 

-  C4g  fi   ! 

K2=I^"   =  la2  x  10   s  (2.25) 

-  2C2g  6   ! 

K4  =  IT   =  3-59  x  10   *  (2.26) 
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It  can  be  seen  that  the  frictional  coefficient  at  the  surface  is 
roughly  three  times  as  large  as  the  internal  coefficient. 

The  mean  and  thermal  vorticity  equations  (2.22)  and  (2.23)  are  the 
prediction  equations  in  the  two  unknowns  q>  and  <JL  .   At  this  point, 
the  geopotential  height  fields  are  restricted  as  follows: 

4h  =   E(y,t)  +  A(y,t)  cos  kx  +  B(y,t)  sin  kx 

+  G(y,t)  cos  mx  +  H(y,t)  sin  mx  (2.27) 

fcj,  ■  F(y,t)  +  C(y,t)  cos  kx  +  D(y,t)  sin  kx 

+  P(y,t)  cos  mx  +  Q(y,t)  sin  mx  (2.28) 

where  E(y,t)  and  F(y,t)  are  the  amplitudes  of  the  zonal  mean  flow; 
A(y,t),  B(y,t),  C(y,t),  and  D(y,t)  are  the  Fourier  amplitudes  of  the 
disturbance  of  wave  number  k,  and  G(y,t),  H(y,t),  P(y,t),  and  Q(y,t) 
are  the  Fourier  amplitudes  of  the  disturbance  of  wave  number  m.   Both 
waves  are  in  the  x-direction  only. 

Since  E  and  F  are  functions  of  y,  barotropic  instability  is  possibl< 
in  this  system  as  well  as  baroclinic  instability.   If  the  definitions 
of  (j)  and  (j)  are  substituted  into  W  =  — —  (k  x  V(j)  ,  we  obtain: 


UH  = 

1  *i 

GM  = 

1     oE 
"    fQoy 

UT  = 

d<l> 

1     °^T 

-  fooF 

GT  = 

1     oF 

'  f    57 

°1 3i 


VM  = 

f    oT" 

o 

_ 

1     oE  . 

VM  = 

f     ox  " 

o 

1        YT 

VT  = 

f    5T" 

_ 

1     ^F   - 

VT  = 

-z     ox 

(2.29) 
(2.30) 

(2.31) 
(2.32) 


where  the  bar  (  )  indicates  the  x-average. 

At  the  north  and  south  boundaries,  the  following  conditions  are 
imposed  initially  and  for  all  time: 
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Uj,  =  11^  =  ^  =  uT  =  0 

at  y  =  0,   y  =  W 

3uM       5uT       ^       SGT        n 

at  y  =  0,   y  ■  W 

where  W  is  the  northern  boundary. 

Substitute  (2.27)  and  (2.28)  into  the  mean  vorticity  (2.22)  and 
thermal  vorticity  (2.23)  equations.   Now  if  we  separate  the  various 
terms,  and  neglect  all  terms  of  wave  number  2k  and  2m  and  combinations 
of  m  and  k,  we  obtain  ten  equations  in  ten  unknowns. 

When  we  equate  the  coefficients  of  the  cos  kx  terms  of  (2.22),  we 
obtain: 

9  ,B2  ,  2,A   k  rS2B  Be  ,  a2D  BF    B3E   nB3£   ,  2,_§E  ,  nB|\, 

^2 
-  j3QBk  -  K4(— 2  -  1c)  (A  -  C)  (2.33) 

The  sin  kx  terms  of    (2.22)   give: 

,B2     ,  2,„       k    r  B2A  Be       b2C  BF    .    .B3E    ,  2>3F    ,    n2,.BE    ,    „Bf, 


°t  *-2  f/   ^„2  By       ^„2  By         ^„3       ^„3 

2 


By  o     By        y        By  By  By" 


+  pQAk  -  K4(^2   -  k2)    (B   -   D)  (2.34) 

By 

The  cos  kx  terms  of    (2.23)    give: 

d    ,h2     .2  ,,2._       k    rBE  B2D    ,   BF  B2B       BE_n  2     2.        dF_„  2     ,2. 
Bt(~2-k  *  )c  =  fQ[B^  ^2  +  B9  ^2  "  By-D(k  ^  >   "  B^B(k  "  ^  > 

3  3  2 

-  £-Id  -  HB^    -   P   Dk  +  K/  (— 9-k2)(A-C) 
oy  ay  oy 

2      B2r 
+  K_    (CkZ  -  2-L)  (2.35) 

2  Sy2 

The  sin  kx  terms  of    (2.23)   give: 

B    ,B2     .2,  ,2.-       kr     BE  B2C       BF  B2A   .   dE_„  2^     2Nj_dFA,.2     ,2. 
^t(~2-k  *  )D  =  fj'  By"  ^2  "  By"  ^2  +  By^  +  ^  >  +  B7A(k  -  ^  } 
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3  3  2  2 

+  ^C  +  ^-|a]    +  (3QCk  +  K4(— 2-k2)(B-D)   +  K2(Dk2-  ^-|)  (2.36) 

dy  dy  dy  dy 

The   terms   independent  of  x  in   (2.22)   give: 

a  *2v-Jl.  ^  t/b    32a  .  „a2p    ^2c, 

d"t  J   2E  "  2f       Zy[\   2   '  *T1       \   2   "   \~2] 
dy  o        J      dy  dy  dy  dy 

m    a  rra2H    „a2o  .  ^fo    „afp, 

+  W   S-y[G7T  -  H7^  +  ^2  -  ^7] 

o  dy  dy  dy  dy 

-h  i  f  -  &  «•»> 

The  cos  mx  terms  of   (2.22)    give: 

d    ,d2  2Xn         m   rd2H  dE    ,    c^Q  dF       ud3E       „d3F        2,TTdE    ,    _d|\ , 

^■m)G=  Va7^    a/5^     a?"  V"m(^    V1 

^2 
-  P  Hm  -  K4(— 2   -  ni  )    (G  -  P)  (2.38) 

dy 

The  sin  mx  terms  of    (2.22)    give: 

d    ,d2  2.u       m   ,     d2G  dE       d2p  dF    .    „d3E    .    J^F        2  ,  3e    ,    pdF,  - 

+  pQGm  -  K   (— 2-  mZ)    (H  -   Q)  (2.39) 

dy 

The  cos  mx  terms  of    (2.23)    give: 

d    ,d2  2     ,2,.       m   rdE  c^Q   ,   dF  d2H       ^E.,   2X  ,  2>      ^F„,   2       ,2. 

3  3  2  2 

-  ^-§Q  -  ^H]    -   P0Qm  +  K4(— 2   -  m2)(G-P)+  K2(Pm2-  §_|)  (2.40) 

dy  dy  dy  dy 

The  sin  mx  terms  of   (2.23)   give: 

d    ,d2        2       2..       m.     dE  d2P       dF  d2G    .   dE^,   2^  ,,2N.   dF_,   2        2, 

d-t(-2-m  "  ^  )Q  =  fo["  57  ^2  "  d7  ^2  +  d7p(m  +  M  )+  d^(m  "  ^  } 

3  3  2  2 

+  ^G  +  ^-|p]    +  PQPm  +  K4(^2-  m2)  (H-Q)+  K^Qm2-  §-§)  (2.41) 

dy  dy  dy  dy 
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The  coefficients   independent  of  x  in   (2.23)    give: 

a    ,d2       ,  .2,„  k       a    r.S2D       32C    ,    22B        32A    . 


5t<-2-  [I  )F  =  2F     §"yt^2  "   BIT  +  CIT  -   °Z2 

2 


dy  o       y     dy  dy  dy  dy' 


~    ^y[GT2  "V2+^2't2+MHP"  GQ)] 
o  dy  dy  dy  dy 

d      ,3e       dFs        ,2  KS2        „     d2F 


+  K4  dV   (Sy"  ■  d7>    "  »     —  "  K2  ^2  <2'42> 

Equations    (2.33)    -    (2.42)   constitute  the  prediction  formulas. 
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III.   PROCEDURE  AND  GOALS 

The  north-south  grid  distance  Ay  was  set  at  200km.  A  time  step, 
At,  of  1/2  hour  was  chosen  to  allow  wind  speeds  up  to  100  m  s  with 
computational  stability. 

Centered  time  differences  were  used  for  all  quantities  except  those 
involving  friction.  The  frictional  terms  were  computed  at  time  t  -  At. 
The  first  time  step  was  a  forward  step.  The  Euler-backward  time  step 
was  introduced  every  24  hours  to  selectively  damp  high  frequency  waves. 
These  waves  are  produced  by  the  separation  of  solutions  at  even  and  odd 
time  steps  caused  by  the  leap-frog  time  scheme. 

The  finite  difference  scheme  used  is  illustrated  below  with  the 
function  F, 

5?  =  2Ay-(Fj+l  -  Fj-1> 

H  =  ~  2  (Fi+l  "  2Fi  +  Fi-1} 
dy2    (Ay)"    J+1      J     J"1   • 

H  =  ~  3[(Fi+2~  2Fi+l+  Fi}-(Fi  ■  2Fi-l  +  Fi-2)] 
dy    2(Ay)J   J+2    J+1   J    J     J  X    J  2 

where  j  is  the  grid  index  and  Ay  the  distance  between  grid  points 

(200km) . 

Equations  (2.33)  through  (2.42)  except  (2.37)  were  solved  by  a 
Gauss  pivotal  elimination  matrix  method  described  by  Richtmyer  (1957). 
Equation  (2.37)  was  solved  by  a  direct  marching  process). 

The  model  was  run  as  a  Beta-plane  experiment  with  periodic  east- 
west  motion  between  north-south  walls  of  separation  W,  4,000  or  8,000 
km. 
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The  8,000  km  experiments  were  likened  to  the  atmosphere.   It  was 
hoped  the  mean  temperature  gradient  from  the  north  to  south  boundary 
at  level  2  would  be  comparable  to  that  of  the  atmosphere  at  500  mb. 
The  source  of  available  potential  energy,  (S*P),  should  yield  a  value 

not  much  less  than  60$  of  the  atmosphere's  total  zonal  generation  of 

-4   -2  -1 
56  x  10  kj  m  s   formalized  by  Dutton  and  Johnson  (1967).   It  should 

be  noted  that  this  model  contains  no  east-west  heating  gradient;  there- 
fore, conversion  of  heat  into  available  potential  energy  (APE)  will  be 
less  than  observed. 

Lastly,  the  winds  at  250  mb  should  exhibit  a  jet  structure  and  the 
750  mb  winds  portray  the  belts  of  zonal  easterlies  and  westerlies. 
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IV.   VERTICAL  MOTION  EQUATIONS 


Vertical  motion  at  500  mb  is  calculated  through  the  thermodynamic 
equation  (2.19). 

COi 


2       ApQ 


rdC       k   /T^E       ^Fv,  .  ^D    ,    k    ,_dE         .^Fn  , 

,s  ^fe  "  f  (D§7  ■  V1     sxn  ^^    f  (C^  "  V] 


4-  cos  mx[^  -  -^  -  H^)]   +  sxn  mx^  4  *  <*§£  -   «g>] 


(4.1) 


all  terms  of  wave  number  2k  or  2m  and  combinations  of  m  and  k  have  been 
neglected. 

The  vertical  motion,  cd,  can  be  written  as 

CD  =  CD  +  CD' 

where  the  average  vertical  motion  is 


Apa 


SF   KS2+  1 

5t  "  ~T  +  2f~  LK  V^  "  **"  T  m  Sy 


-  [k  x.(BC  -  AD)  +  m  ^„(PH  -  GQ) ]  }      (4.2) 
o 


again,  the  barred  quantity  means  x-averaged. 

The  disturbance  vertical  motion,  cd',  is  given  by 

cd'    =  cd1    (cos  kx)   +  cd'  (sin  kx)   +  cd'  (cos  mx)    +  cd1  (sin  mx) 
c  s  c  s 

where 

cd^cos  kx)    =  ^     (cos  kx[^  -  ^-(D^  -  B^)]  J 

i  /    •      in-     2        (     .      ,     r^D    ,    k    ._dE        3FS1 

a>;("ia  kx)  =  ^  ( sin  kxfe +  r(%  -  v1 


1/  N     -       2 

cd   (cos  mx)    =  — t— 
c  aAp 


rdP       m   „3E       TTdFN, 

COS      ^[-.—(Q^      .      1^)] 

f  _ .  r^O    .    m    /T,dE        „dF. 


es(sin  ra)  =  Hf  (sin  »f5t +  ^(p§7  -  V1 


(4.3) 

(4.4) 

(4.5) 

(4.6) 
(4.7) 
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V,   COMPUTATION  OF  AVERAGE  NORTH- SOUTH  WIND 

The  average  north-south  flow  computations  coupled  with  the  average 
vertical  motion  provides  one  with  a  fast  method  of  calculating  the  cir- 
culation in  the  meridional  plane;  i.e.  Hadley  or  Ferrel  cells. 

If  we  apply  the  continuity  equation  at  level  1,  we  obtain 
& 


V  •  \V  (5.1) 


^Vl      ^°2 
by  dp 

Thus,  the  finite  difference  form  is 


where   (v..).    =0   . 


(5.2) 


-  /W,         faJi     +     ^o)  4-1.1 

<Vl  ■   «Vj   -  I   '     2    1   2        2   1+1'  <5-3> 
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VI.   ENERGY  AND  ENERGY  TRANSFORMATION  EQUATIONS 

The  energy  equations  were  obtained  from  Phillips  (1956). 
Mean  Potential  Energy  per  unit  mass, 

P  sJL- o  |(1j  "  ?o)2dy  (6.1) 

8Wf  *  J      L         J 
o 

Disturbance  Potential  Energy  per  unit  mass, 


{<*;-*; 


,2    (        , 

,       ro)  dy  (6.2) 

8Wf  Z  ' 
o 


Mean  Kinetic  Energy  per  unit  mass. 


1  (\-2  .  -5 

^-J(ui+U: 


"  4F  I  (U1  +  U3)dy  (6'3) 


Disturbance  Kinetic  Energy  per  unit  mass. 


f[(v^: 


^-o  I  [(V^)Z  +  (vk)Z]  dy  (6.4) 

4Wf 
o 


If  we  substitute  the   (J)     and    <JL  into   the  above  energy  equations, 


we  obtain 

,2  N 


P  =  -*=—,     2       F2 

2Nf  j=0     J 

o        J 

2  N  2  N 

P'(k)    =  -^—  2      (C2  +  D2)  P'(m)    =  -^— ?     2        (P2  +  Q2) 

4Nf  1=0     J  J  4Nf  j=0        J  J 

o  J  o        J 

k  -  _i v      rr^2  +  A2i 

K"2Nfo2    h^  +  W 

o         J 
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■■«^L^+^  +  *  +  *5  +  *' 


where  the  subscripts  k  and  m  stand  for  wave  number,  and  N  is  the 
number  of  grid  increments  such  that 

W  =  NAy 
Following  the  procedure  of  Phillips  (1956),  the  energy  flow  diagram 
of  Figure  2  should  make  the  physical  interpretation  of  the  energy  trans- 
formation equations  clear. 


Cs-p)- 


(P-K) 


_J_ 


T 


K    —fc-t,) 


(P'P') 

1 


(K'.R) 


Cp'.k') 


-<*'*<) 


Fig.  2.   Energy  Flow  Diagram 
The  flow  of  energy  is  in  the  direction  of  the  arrows  if  the  asso- 
ciated transformation  (  •  )  is  positive. 

The  energy  transformation  equations  are: 

Ll2H  N 

(S-P)    =  -*-£=-     S        (S.    •   F.) 

2Nf  j=0        J  J 


(P'P')(k) 


.jA 


,dF, 


^-     2        [(A-D  -   B-C)      •    (^-)    ] 

2Nf  3       j=0  j        ^  J 

o  J 

N 


(P'-K')(k)    = 


2Ap  N 


£        [C.cd'  .(cos  kx)   +  D.to'  .(sin  kx)  ] 
j=0      J  CJ  J  SJ 
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+  £L\       (pd  C  _  A^2D  +  p^2A  _  C^2B.  , 
y  -*   dy     dy     dy     dy   -" 

o  J 

K2    N     2  2    2       B2C    S2D 
(K'-O(k)  =-^-?  2   [k2(C2  +  D2).  -  <(£-£  +  d2-|)  ] 

2       2NfQ2   j=0  J     byZ  ay2  J 

-K,    N  -.2^   ^2,, 

(K.V=  -4l    (E-F)   •  ^-^f). 
4     N  fQ2   j=0        J    -Sy2   Sy2  J 

K    N   (     ^2    >2 

(K'-T,)(k)  =  V  2   -A .[(^-f  -  ^).  -  k2(A-C).] 

4       2NfQ2   j=o(   J   Sy2   ^2  J  J 

-B.[(^|-^f)i  -k2(B-D)J 
J   ^y2   By2  J  J 

The  energy  transformations  of  the  m  wave  disturbance  have  been 

omitted  for  brevity.   By  replacing  k,  A,  B,  C,  and  D  with  m,  G,  H,  P, 

and  Q,  the  entire  transformation  set  will  be  complete. 
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VII.   INITIAL  CONDITIONS 

A  natural  initial  state  would  be  one  of  no  motion  and  constant  tem- 
perature.  Then  the  system  could  run  until  a  steady  state  is  reached  or 
some  type  of  instability  occurs.   However,  to  save  computer  time,  the 
initial  state  was  chosen  to  be  near  a  critical  state  for  development  of 
the  most  unstable  wave. 

When  the  mean  flow  model  was  allowed  to  run  with  no  disturbance, 
the  wind  field  reached  a  steady  state;  the  Reynolds  stresses  balanced 
the  heat  input.   This  steady  state  was  comprised  of  a  meridional  Hadley 
cell  superimposed  on  the  zonal  flow.   This  Hadley  cell  mean  flow  was 
utilized  as  the  initial  condition. 


If 


we  set  «r-  =  0  and  A-D,  G-Q  =  0 


5  <^E   ^F 
the  F  equation  yields  0  =  K4  ^y(^  -  ^) , 

2 

the  E  equation  yields  0  =  -|J2Kh(^  -  y)/W  -  K.  2JL  . 

2     .   2  ^y2 

The  solution  for  F  is: 

2      2    3 

*    K2   ^4    6W;  U#i; 

This  solution  requires  that  E  =  F  which  implies  that  (JL  =  (j)  .   In  this 
case 

u3  =  0  and  U;L  =  2uT  =  2^  . 

Pedlosky  (1964)  derived  a  simple  stability  criterion  for  the  mean 
flow.   If  the  potential  vorticity  gradient  («A  changes  sign  in  either 
the  y-direction  or  the  vertical  instability  is  possible. 
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For  this  model  dq/3y  is  given  by 
v2- 


'      -       °  ui  ,  u2 


s  3       3%   .,2  _    _ 

5^  ■  p0  -  ^r  -  £-<»!  -  V  <7-3> 

Although  u  =  0  initially,  easterlies  and  westerlies  develop  with 

time  as  observed  in  the  earth's  atmosphere.   However,  because  u  <<  u  , 

we  shall  neglect  u_  for  this  discussion. 

If  we  neglect  u  in  equation  (7.3),  we  see  that 
3  2 

is  the  critical  value  of  the  mean  flow.  Above  that  value,  instability 

may  occur  for  some  wave  number.   If  we  solve  for  u-    .  .   ,  a  value  of 

h   .  .   ,  may  be  derived  through  the  Hadley  cell  initial  condition  (7.1). 

In  this  paper 

h      .„.      .    =  0.0011 
critical 

The  disturbance  initial  conditions   are 

A  =  200  sin  TTy/w 

B,    C,    D     =0 

G  =  200  sin  TTy/w  after  100  days 

H,    P,   Q     =0 
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VIII.   RESULTS 

A.  HIERARCHY  OF  UNSTABLE  WAVE  NUMBERS 

The  objective  of  the  first  experiment  was  to  find  a  hierarchy  of 

wave  numbers  in  order  of  wave  growth  for  a  particular  mean  flow.   The 

mean  flow  was  given  by  the  Hadley  cell  solution  with  h  =  0.004  kj  ton 

s   and  it  was  held  constant  with  time.   Since  the  mean  flow  did  not 

vary  with  time,  the  equations  for  the  disturbance  became  linear.   Each 

case  was  run  until  expotential  growth  was  achieved.   Figure  3  shows  the 

graphs  of  wave  number  vs  wave  growth  rate  for  two  sets  of  conditions: 

B  =  0  and  P  =  1.67  x  lO"11,  and  W  =  4,000  km  and  W  =  8,000  km.  The 
ro         o  *  ' 

growth  rate  shown  is 

1     .  Amp(t2)      ? 

growth  rate  =  ^-  In  t^^l  x  10 

where  Amp(t)  =  amplitude  of  the  disturbance  at  time  t. 

Pedlosky  (1964)  has  computed  the  growth  rate  for  this  wind  profile 
using  different  values  for  various  parameters. 

Figure  3  exhibits  how  the  Beta  term  stabilizes  the  longer  waves. 

B.  INTERACTIONS  BETWEEN  WAVE  AND  MEAN  FLOW 

The  second  set  of  experiments  were  devised  to  test  the  growth  of 
the  most  unstable  amplifying  wave,  allowing  interactions  with  the  mean 

flow,  as  a  function  of  heating.   All  experiments  were  conducted  using 

-4      -1  -1 
Palmen's  exchange  coefficient  of  225  x  10   ton  m   s   and 

heat  =  h   _  .   -  (0.001  kj  ton   s   )  as  the  initial  condition  for  the 
critical 

Hadley  cell. 
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Each  experiment  was  run  for  300  days.   The  k  wave  was  introduced 
immediately;  the  m  wave  (if  existed)  was  introduced  at  100  days. 

The  following  tables  were  set  up  to  briefly  describe  the  character- 
istics of  each  wave  number  by  way  of  symbolic  nomenclature.   The  sym- 
bols are: 

D  =  wave  dies  immediately 

N  =  a  non-steady  wave,  disturbance  energy 

oscillates  up  and  down  as  the  wave  interacts 
with  the  mean  flow. 

SW  =  a  steady  wave,  total  disturbance  energy 

2  -2 

is  =  (   )  m  s 

Peak  -  the  peak  amount  of  disturbance 

2  -2 
energy  =  (   )  m  s   at  day. 

E1 (k)  vs  E'(m)  -  comparison  of  k  and  m  wave 

disturbance  energy  over  the  200  day  period 

when  both  waves  are  introduced. 
E'(m)  >  E'(k)   -  if  che  m  wave  disturbance  energy 

did  exceed  that  of  the  k  wave,  then  the  m  wave 

2  -2 

disturbance  energy,  E'(m)  =  (   )  m  s 

on  day. 

-2-1 

S»P  =  (   )  kj  m  s   is  the  average  amount  of  conversion  of  heat 

to  APE  usually  computed  after  the  mean  flow  has  come  to  some 

degree  of  equilibrium  with  the  disturbance.   If  the  system  has 

not  run  long  enough  to  reach  the  equilibrium,  then  the  value  is 

given  which  best  represents  the  average  S*P  thus  far  incurred. 
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C.   U-  AND  U_  VERSUS  HEATING  FOR  THE  STEADY  WAVE  EXPERIMENTS 

The  table  below  shows  the  dependence  of  the  zonal  flow  upon  the 
net  heating  rate  for  the  Steady  Wave  experiments. 


heat   (kj 

ton 

•1 

-1 
s 

exp# 

uT(m 

s"1) 

u^m  s      ) 

0.004 

3 

6.9 

14.1 

0.008 

4 

7.1 

15.1 

0.012 

5 

7.2 

15.6 
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Table  5.11  lists  the  steady  state  values  of  the  zonal  flow.   The 

rather  small  increases  in  u  with  heat  is  best  explained  by  the  rapid 

i. 
baroclinic  response  to  the  wind  sheer  decaying  the  unstable  zonal  flow. 


A  BAROCLINIC,  STEADY  WAVE 


=  1.67  x  10"11 
=  4,000  km 
Experiment  #5  \   k    =7  /  ,  exhibits  a 

heat  =  0.012 


very  interesting  feature  in  this  model:  a  steady,  baroclinic  wave. 
Figure  4  is  a  graph  of  mean  kinetic  energy,  K,  mean  potential  energy, 
P,  mean  energy,  E,  and  total  energy,  E,  for  experiment  #5.   Initially, 
the  mean  flow  grows  rapidly,  as  shown  on  Figure  4,  with  negligible 
disturbance  energy.   At  the  40th  day,  the  zonal  mean  flow  has  reached 
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its  peak  and  decreases  rapidly  while  the  disturbance  begins  to  grow. 
By  the  115th  day,  the  entire  atmosphere  has  reached  a  steady  state. 

Initially,  the  maximum  u-  =  4.9  m  s   and  u  =  2.4  m  s   .   At  the 
20th  day,  the  maximum  u  =  14.2  m  s   (henceforth,  whenever  zonal  mean 
flow  is  mentioned,  it  will  mean  the  maximum  speed  in  the  north-south 
column)  and  u  =  7.5  m  s   .   The  40th  day  exhibited  the  maximum  mean 
flow  energy,  due  to  the  maximum  P,  with  u-  =  21.7  m  s   and  u  =  11.1 
m  s   .   However,  the  winds  continued  to  increase  until  the  46th  day, 
as  attested  to  by  the  K  graph,  with  u_  =  25.3  m  s   and  u  =  10.9  m  s 

Figure  5  (E'(k))  portrays  how  the  wave  grew  beginning  at  the  32nd 
day.   By  the  38th  day,  the  wave  responds  to  the  strong  baroclinic  in- 
stability, as  shown  by  the  large  u  ,  and  grows  expotentially.   Baro- 
tropic  damping  (K' »K)  occurs  exclusive  of  any  barotropic  growth  the 
entire  cycle  of  the  wave,  and  exhibits  large  oscillations  during  the 
rapid  growth  of  the  wave. 

The  wave  continued  to  grow  by  baroclinic  instability  until  the  56th 

day  when  the  wave  reached  its  peak  energy  with  a  disturbance  energy,  E', 

2  -2 

=  123  m  s   .   By  this  time,  the  thermal  wind  had  decreased  from  a  peak 

of  11.2  m  s   on  the  42nd  day  to  7.5  m  s   on  the  56th  day,  destroying 
the  source  of  the  baroclinic  growth. 

On  the  38th  day,  when  the  wave  began  its  large  baroclinic  growth, 
barotropic  damping  begins  to  reach  its  greatest  value  as  seen  in  Fig.  5. 
Initially,  as  the  barotropic  damping   oscillates  about  its  mean  of 
1.4x  10  kj  ton  s, the  baroclinic  instability  is  by  far  dominant.   How- 
ever, by  the  58th  day,  barotropic  damping  plus  the  loss  of  eddy  kinetic 
energy,  K1 ,  by  frictional  dissipation  becomes  larger  than  PI#K';  hence, 
baroclinic  effects  must  be  diminishing. 
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As  the  barotropic  damping  reaches  steady  state,  it  appears  to  come 
into  equilibrium  with  baroclinic  growth,  hence,  the  steady  state  of  the 
disturbance  energy. 

Figure  6  is  a  graph  of  the  conversion  of  net  heating  to  available 

potential  energy  (S*P).   The  graph  shows  the  rapid  rise  of  S#P  to  a 

-4    -2  -1  -4      -1 

peak  of  16  x  10  kj  m   s   in  a  vertical  column  (1.6  x  10  kj  ton 

s   )  at  40  clays,  then  its  decrease  to  the  steady  state  value  of  8.8 

-4    -2  -1 
x  10  kj  m   s   by  the  104th  day. 

At  the  end  of  the  300th  day  run,  the  entire  atmosphere  has  been  at 
a  steady  state  equilibrium  since  the  115th  day.   The  energy  conversion 
terms  feeding  P,  P1,  K,  and  K1  were  constant  with  respect  to  time,  as 
were  the  energy  terms;  therefore,  the  conversion  terms  should  balance 
for  each  energy  term.   Figure  7  is  a  pictorial  display  of  the  300th  day 
energy  flow  after  Phillips  (1956),   the  only  conversion  not  balancing 
being  the  (K* T  ) . 

In  summary,  the  steady  wave  is  achieved  by  the  conversion  of  P  to 
K'  which  is  ultimately  balanced  by  frictional  losses  and  barotropic 
damping  through  the  action  of  the  Reynolds  stresses. 

E.   A  N0N- STEADY  WAVE  WITH  BAROTROPIC  EFFECTS 

0  =  1.67  x  10"11 
w°  =  8,000  km 
Experiment  #12l   k  =  6  /  ,  exhibits  another 

\heat  =  0.016 

quite  interesting  feature,  a  wave  initially  amplified  by  baroclinicity 

which  then  grew  by  barotropic  effects. 

Figure  8  is  a  graph  of  the  mean  flow  energies  and  total  energy  for 

experiment  #12.   Notice  the  large  and  small  oscillations  in  each  energy 

graph  and  how  the  total  energy  graph  is  smoother  than  the  mean  flow 
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energy  graph.  This  is  due  to  the  disturbance  drawing  energy  from  the 
mean  flow  and  vice  versa.   So  when  the  disturbance  decreases,  the  mean 
flow  increases  and  vice  versa  which  tends  to  smooth  out  the  total 
energy  graph. 

On  the  20th  day,  the  wave  grew  rapidly  by  baroclinic  effects.   At 
that  time,  u  =  9.3  m  s   and  u_  =  18.1  m  s   while  on  the  30th  day, 
uT  =  10.9  m  s   and  u  =  30.0  m  s   .   Since  only  barotropic  damping  was 
occurring  and  the  mean  flow  increased  since  the  initial  10  days  of 

growth,  baroclinic  effects  amplified  the  wave  to  its  initial  peak  of 

2  -2 
69.4  m  s   at  30  days. 

Figure  9  is  a  graph  of  the  disturbance  energy  and  the  barotropic 
instability.   It  may  be  seen  that  the  barotropic  growth  term  feeds  or 
damps  the  disturbance.   At  day  30,  the  wave  has  peaked  due  to  baro- 
clinic growth.  From  the  31st  to  the  33rd  day,  the  baroclinic  instability 
subsided  and  the  barotropic  damping  destroyed  over  half  of  the  wave 

energy  and  increased  the  mean  flow  to  u  =  15.3  m  s   and  u  =  38.2  m 

-1 
s 

The  34th  is  the  first  day  that  barotropic  instability  occurred  and 

2 

although  it  only  lasted  for  two  days,  the  disturbance  peaked  at  77  m 

-2 

s   .   On  the  37th  day,  barotropic  damping  again  diminished  the  disturb- 
ance until  the  40th  day.  On  the  37th  day,  Up  =  14.1  m  s   and  u..  = 
37.9  m  s   while  on  the  40th  day,  u  =  17.2  m  s   and  u  =  42.6  m  s 

To  further  supplement  this  discussion,  the  K  graph  is  given  on  the 
bottom  of  Figure  8,  without  its  energy  scale.   The  graph  helps  to  por- 
tray the  barotropic  damping  effect,  gain  of  mean  flow  energy  to  the 
disturbance  and  vice  versa. 
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On  the  40th  day,  the  barotropic  instability  causes  explosive  ampli- 
fication of  the  wave.   However,  in  the  42nd  and  43rd  days,  barotropic 
damping  occurs  and  explosive  baroclinic  growth  of  approximately  five 
times  the  average  P' *K'  takes  over. 

This  baroclinic  growth  can  be  seen  on  the  small  growth  of  the  K 
graph  from  the  initial  loss  due  to  barotropic  growth  from  the  40th  to 
the  42nd  day  to  the  small  gain  from  the  42nd  to  the  44th  day. 

From  the  44th  to  the  53rd  day,  growth  is  a  combination  of  baro- 
tropic and  baroclinic  effects.   Possibly  the  barotropic  term  is  domin- 
ant, as  with  each  barotropic  damping,  the  disturbance  graph  changes 
slope  or  declines,  hence  a  greater  dependence  on  barotropic  growth. 
The  mean  flow  has  declined  greatly  due  to  both  a  baroclinic  and  baro- 
tropic drain  on  it,  as  shown  by  the  44th  day  u  =  14.2  m  s   and 
u1  =  40.5  m  s   while  on  the  53rd  day,  u  =  11.9  m  s  and  u  =  23.9  m  s   . 

From  the  53rd  to  the  82nd  day,  the  wave  damped  for  the  most  part 
due  to  barotropic  effects,  although  a  small  amount  of  baroclinic  damp- 
ing does  occur. 

The  disturbance  peaks  on  the  56th  day,  due  to  barotropic  and  baro- 
clinic growth,  with  the  K  graph  reaching  its  minimum  on  the  55th  day. 
Again,  baroclinic  effects  do  create  growth  in  the  wave  as  seen  when 
barotropic  damping  is  occurring  and  the  disturbance  pauses  from  its 
decay  to  grow  for  a  day,  then  continue  its  decay. 

Figure  10  is  included  to  observe  the  variation  of  S*P  term  in  a 
non-steady  wave.  Figure  11  is  included*  to  show  the  many  oscillations 
of  the  P#K  term  (thin  line).   An  interesting  phenomenon  noted  on  Figure 
11  is  the  comparison  of  the  P#K  conversion  with  the  K  term  (heavy  line, 
without  energy  scale) .   The  two  terms  tend  to  counteract  each  other 
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rather  than  reinforce  each  other.   In  the  same  fashion,  K' *K  also 
counteracts  the  P»K  conversion  with   K,#K  on  the  average  dominating 
P»K.   Hence  with  barotropic  damping,  (K' *K  >  0)  the  P#K  conversion 
usually  becomes  negative.   The  reverse  also  occurs. 

Figure  12  is  a  graph  of  the  u.  flow  with  respect  to  time.   The 
graph  begins  when  the  disturbance  and  mean  flow  reach  a  state  of 
equilibrium  after  the  initial  growth  period  of  the  disturbance  on  the 
mean  flow.   The  graph  runs  from  80  to  110  days.   The  easterlies  present 
at  250  mb  are  attributed  to  the  fact  that  no  horizontal  diffusion  is 
present  in  the  vorticity  equation. 

F.   A  TWO  WAVE  EXPERIMENT 

0   =  1.67  x  10"11 

.  w°  =  8,000  km 
Experiment  #12 (  ,   _  ,  ,  is  one  of  the 

m  -  7 
heat  =  0,012 

few  experiments  conducted  where  E'(m)  >  E'(k).   The  explanation  for  the 
dominance  of  the  m  wave  over  the  k  wave,  in  opposition  to  the  results 
of  section  V.  A.,  is  attributable  to  the  modification  of  the  mean  zonal 
flow  by  the  k  wave.   Figure  13  is  a  graph  of  E'(k)  and  E' (m)  versus 
time.   Figure  14  is  a  graph  of  the  K1 *K  conversions  for  each  wave.   It 
should  be  noted  that,  on  the  average,  both  conversions  act  in  the  same 
direction  in  quite  varying  degrees. 

The  Pedlosky  criteria  (7.2  and  7.3)  proved  to  be  a  necessary  condi- 
tion for  baroclinic  and  barotropic  instability.   The  baroclinic  criteria 
was  met  for  the  initial  condition  and  throughout  the  forecast.   The 
barotropic  criteria  was  met  continuously  after  the  34th  day.   At  no  time 
did  barotropic  growth  occur  when  the  Pedlosky  criteria  was  not  met. 
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Between  the  162nd  and  164th  day,  the  E'(m)  becomes  larger  than  the 
E'(k)  and  remains  such  throughout  the  run.   At  this  time,  both  waves 
are  experiencing  barotropic  instability.   However,  the  m  wave  is  grow- 
ing and  the  k  wave  is  decaying. 

The  time  frame,  from  162  to  164  days,  is  explored  further  on 
Figures  15,  16  and  17.   Figure  15  and  16  is  a  graph  of  P' »K*  (bold  line) 
and  K1 *K(thin  line)  for  waves  k  and  m  respectively.   A  unique  differ- 
ence is  apparent  between  the  two  waves.   The  baroclinic  conversion 
P'^K1  never  approaches  zero  with  the  m  wave  but  does  many  times  with 
the  k  wave,  and  even  becomes  negative.   When  the  baroclinic  term  be- 
comes negative,  baroclinic  damping  draws  from  K'.   Usually,  the  P'*K' 
conversion  was  positive,  thus  the  waves  were  experiencing  baroclinic 
growth.   On  day  164,  both  waves  are  experiencing  barotropic  growth; 
however,  the  k  wave  is  undergoing  baroclinic  decay  and  the  m  wave  is 
experiencing  baroclinic  growth.   The  energy  flow  chart  is  shown  below 
to  help  explain  the  growth  of  the  m  wave  and  decay  of  the  k  wave. 


[I 


1.41 


088 

r  I 

LpJ^0.16-4_KJ-*0.38  1    P'   [-o,80-4~KH-->  0 


44 


026 


032 


net    change    to    K'(k)  net    change    to    Kf(m) 

-1        -1  -1      -1 

=    +.61kj    ton         s  =    +.92kj    ton         s 

Possibly  this  difference   in  net   inflow  to   the  K1    term  for  each  wave 

gives   the  growth  advantage  to   the  m  wave  for  a  ten  day  period. 
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If  we  look  at  the  waves  on  the  164th  day,  we  find  an  interesting 
feature  of  the  barotropic  growth.   Although  the  k  wave  K' *K  is  largest 
(negative),  the  m  wave  is  experiencing  the  growth  through  baroclinic 
instability.   Since  the  conversion  terms  are  the  meridional  average, 
it  is  possible  the  barotropic  growth  of  the  m  wave  in  a  sector  is  sig- 
nificantly greater  than  barotropic  growth  of  the  k  wave  in  a  similar 
region.   This  possibility  is  explored  in  Figure  17,  a  graph  of  the  ft. 
phase  angle  of  the  k  and  m  waves  and  u  .   The  wave  is  experiencing 
barotropic  growth  if  the  tilt  of  the  wave  is  opposite  the  shear  of  the 
mean  flow.   Figure  17  depicts  the  k  wave  growing  barotropically  south 
of  the  zonal  maximum  of  u  to  the  easterly  jet,  and  the  m  wave  growing 
barotropically  north  of  the  zonal  maximum  to  about  y  =  6,000  km,  thus 
demonstrating  the  greater  barotropic  instability  of  the  k  wave. 
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IX.   CONCLUSIONS 

In  order  to  facilitate  physical  interpretation  of  the  results,  the 
model  was  restricted  to  a  mean  flow  and  one  or  two  waves.   A  steady 
state  solution  which  was  independent  of  x  was  obtained.   This  flow  was 
of  the  Hadley  type  with  zero  mean  velocity  at  level  3  and  a  parabolic 
wind  profile  at  level  1.   The  stability  of  this  mean  flow  was  investi- 
gated by  performing  initial  value  integrations  with  the  linearized 
frictionless  equations.   The  calculations  were  continued  until  expon- 
ential wave  growth  was  achieved.   At  this  point,  the  growth  rates  for 
a  range  of  wave  numbers  were  calculated.   The  growth  rate  of  the  waves 
was  determined  as  a  function  of  wall  separation  and  |3.   The  growth  rate 
was  reduced  by  the  (3  effect  and  increased  greatly  by  large  wall  separa- 
tions. 

The  nonlinear  frictional  equations  were  integrated  from  an  initial 
state  which  consisted  of  the  Hadley  mean  flow  upon  which  was  super- 
imposed the  wave  number  of  maximum  growth  rate.   For  each  combination  of 
heating,  P,  and  wall  separation,  a  300  day  prediction  was  made  with  the 
optimum  wave  number  from  the  linear  study  previously  mentioned.   Each 
of  these  forecasts  was  repeated  but  at  100  days  a  new  wave  was  intro- 
duced.  The  purpose  here  was  to  determine  the  largest  amplitude  wave 
which  would  exist  on  the  mean  flow,  the  latter  having  been  modified  by 
nonlinear  interactions  by  the  first  wave.   A  total  of  45  experiments 
were  made  as  a  function  of  wall  separation,  heating,  |3,  and  wave  number. 

Constant  amplitude  propagating  disturbances  were  obtained  in  every 
case  where  the  wall  separation  was  4,000  km.   These  experiments 
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contained  cases  where  (3  =  0  and  M  0  as  well  as  a  large  range  in  heat- 
ing rates.   The  steady  waves  gained  energy  through  baroclinic  conver- 
sion processes  and  lost  energy  through  barotropic  decay  and  friction. 

The  remaining  experiments  with  the  wall  separation  of  8,000  km  pro- 
duced non-steady  waves  with  fluctuating  baroclinic  and  barotropic  in- 
teractions.  The  roles  of  the  different  types  of  instability  in  the 
growth  and  decay  of  the  waves  were  easily  distinguished.   When  one  in- 
stability was  causing  strong  growth,  the  other  instability  was  usually 
a  minimum  growth  conversion  or  was  decaying  the  disturbance.   The  baro- 
tropic and  baroclinic  conversions  tended  to  be  out  of  phase  with  one 
another,  one  a  maximum  -  the  other  a  minimum.   The  P*K  conversion, 
however,  tended  to  be  a  constraint  on  the  K  term  and  adjusted  to  the 
barotropic  terra.   With  strong  barotropic  damping,  the  P*K  was  negative, 
while  with  strong  barotropic  growth,  the  P*K  conversion  was  positive. 
Hence  with  strong  barotropic  decay,  the  Ferrel  cell  became  dominant  over 
the  Hadley  cell. 

Over  the  long  term  mean,  the  waves  obtained  a  net  gain  of  energy 
through  P,#K'  and  a  net  loss  of  energy  through  K' #K.   The  mean  tempera- 
ture gradient  from  the  south  wall  to  the  north  wall  and  the  S#P  conver- 
sion were  directly  related  to  the  heating  and  wall  separation.   With 
the  largest  heating  and  wall  separation,  the  temperature  gradients  were 
somewhat  greater  than  the  north-south  500  mb  temperature  gradients  in 
the  earth's  atmosphere.   This  is  related  to  the  larger  than  normal 
heating  used  to  offset  the  lack  of  east-west  heating  gradient.   Even 
with  the  largest  heating  and  wall  separation,  the  conversion  of  heat  to 

available  potential  energy,  S»P,  is  only  about  65$  of  the  average  for 

-4    -2  -1 
our  atmosphere  of  56  x  10  kj  m   s   as  estimated  by  Dutton  and 
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Johnson  (1967).   This  underestimate  is  also  the  result  of  the  lack  of 
east-west  heating. 

The  zonal  mean  flow  at  level  1  exhibited  a  definite  westerly  jet 
in  the  mid- latitudes,  which  sometimes  split  into  a  double  jet.   Due  to 
the  lack  of  a  horizontal  diffusion  term,  the  horizontal  gradients  were 
sometimes  quite  large  and  occasionally  an  easterly  jet  would  appear  for 
a  short  time.   The  zonal  mean  flow  at  level  3  exhibited  the  belts  of 
zonal  easterlies  and  westerlies  typical  of  the  earth's  atmosphere. 

In  this  study,  we  found  with  a  wall  separation  of  4,000  km  only 
steady  baroclinic  waves,  and  with  a  wall  separation  of  8,000  km 
fluctuating  baroclinic  and  barotropic  interactions  yield  non-steady 
waves.   It  was  determined  that  the  predominate  wave  number  in  the  fully 
developed  flow  was  in  every  case  close  to  the  value  determined  by  the 
linear  stability  analysis.   In  the  actual  atmosphere,  where  walls  are 
not  present,  barotropic  instability  may  be  important  in  giving  day  to 
day  fluctuations  in  the  mean  flow. 


63 


LIST  OF  REFERENCES 


1.  Dutton,  J.  A.  and  D.  R.  Johnson,  1967:   The  theory  of  available 

potential  energy  and  a  variation  approach  to  atmospheric 
energetics.   Advances  in  Geophysics,  Vol.  12,  333-436. 

2.  Fultz,  D.  et  al,  1959:   Studies  of  thermal  convection  in  a  rotating 

cylinder  with  some  implications  for  large  scale  atmospheric 
motions.  Meteorological  Monographs,  Boston,  American  Meteoro- 
logical Society,  4,  21,  104  pp. 

3.  Herring,  J.  R. ,  1963:   Investigation  of  problems  in  thermal  convec- 

tion.  Journal  of  the  Atmospheric  Sciences,  Vol.  20,  325-338. 

4.  ,  1964:   Investigation  of  problems  in  thermal  convec- 


tion: Rigid  boundaries.   Journal  of  the  Atmospheric  Sciences, 
Vol.  21,  277-290. 

5.  Hide,  R. ,  1953:   Some  experiments  on  thermal  convection  in  a  rotat- 

.  ing  liquid.   Quarterly  Journal  of  the  Royal  Meteorological 
Society,  79,  161. 

6.  Kasahara,  A.  and  W.  M.  Washington,  1967:   NCAR  global  general  cir- 

culation model  of  the  atmosphere.   Monthly  Weather  Review,  Vol. 
95,  No.  7,  389-402. 

7.  ,  1971:   General  circulation  experiments  with  a  six- 
layer  NCAR  model,  including  orography,  cloudiness,  and  surface 
temperature  calculations.   Journal  of  the  Atmospheric  Sciences, 
Vol.  28,  No.  5,  657-701. 

8.  Manabe,  Syukuro,  1969:   The  atmospheric  circulation  and  the  hydro- 

logy of  the  earth's  surface.  Monthly  Weather  Review,  Vol.  97, 
No.  11,  739-774. 

9.  Manabe,  Syukuro  and  J.  Smagorinsky,  1967:   Simulated  climatology  of 

a  general  circulation  model  with  a  hydrologic  cycle:  II. 
Analysis  of  the  tropical  atmosphere.   Monthly  Weather  Review, 
Vol.  95,  No.  4,  155-169. 

10.  Matsuno,  J.,  1966:   Numerical  integrations  of  the  primitive  equa- 

tions by  a  simulated  backward  difference  method.   Journal  of 
Meteorological  Society  of  Japan,  44,  76-84. 

11.  Mintz,  Y.,  1965:   Very  long-term  global  integration  of  the  primitive 

equations  of  atmospheric  motion.   WMO  Technical  Note  No.  66, 
"WMO-IUGG  Symposium  on  Research  and  Development  Aspects  of  Long 
Range  Forecasting,  Boulder,  Colo.,  1964,"  Geneva,  1965,  141-167. 


64 


12.  Palmen,  E.,  1955:   On  the  mean  meridional  circulation  in  low  lati- 

tudes of  the  northern  hemisphere  in  winter  and  the  associated 
meridional  and  vertical  flux  of  angular  momentum.   Societas 
Scientiarum  Fennica,  Commentations  Physica  Mathematicae,  Vol. 
17,  No.  8,  1-33. 

13.  Pedlosky,  J.,  1964:   The  stability  of  currents  in  the  atmosphere 

and  the  ocean.   Journal  of  Atmospheric  Science,  Vol.  21,  201- 
219,  342-353. 

14.  Phillips,  N.  A.,  1956:   The  general  circulation  of  the  atmosphere. 

A  numerical  experiment.   Quarterly  Journal  of  the  Royal 
Meteorological  Society,  Vol.  82,  No.  352,  123-164. 

15.  Richtmyer,  R.  D. ,  1957:   Difference  Methods  for  Initial  Value  Prob- 

lems.  Interscience  Publishers,  Inc.,  New  York. 

16.  Riehl,  H. ,  et  al,  1951:   The  North-East  trade  of  the  Pacific  Ocean. 

Quarterly  Journal  of  the  Royal  Meteorological  Society,  Vol.  77, 
No.  334,  598-626. 

17.  Rossby,  C.  G. ,  and  R.  B.  Montgomery,  1935:   The  layer  of  frictional 

influence  in  wind  and  ocean  currents.   Papers  in  Physical 
Oceanography  and  Meteorology,  Vol.  3,  No.  3,  101. 

18.  Smagorinsky,  J.,  1963:   General  circulation  experiments  with  the 

primitive  equation,  I.   The  basic  experiment.  Monthly  Weather 
Review,  Vol.  91,  No.  3,  99-165. 

19.  Smagorinsky,  J.,  S.  Manabe  and  J.  L.  Holloway,  Jr.,  1965:   Numerical 

results  from  a  nine-level  general  circulation  model  of  the 
atmosphere.   Monthly  Weather  Review,  Vol.  93,  No.  12,  727-768. 


65 


INITIAL  DISTRIBUTION  LIST 


No.  Copi< 


1.  Defense  Documentation  Center  2 
Cameron  Station 

Alexandria,  Virginia  22314 

2.  Library,  Code  0212  2 
Naval  Postgraduate  School 

Monterey,  California  93940 

3.  Associate  Professor  Roger  T.  Williams,  Code  51  10 
Department  of  Meteorology 

Naval  Postgraduate  School 
Monterey,  California  93940 

4.  Lieutenant  Glenn  Curtis  Trumbower,  USN  5 
Fleet  Numerical  Weather  Central,  Naval  Postgraduate  School 
Monterey,  California  93940 

5.  Naval  Weather  Service  Command  1 
Washington  Navy  Yard 

Washington,  D.  C.   20390 

6.  Officer  in  Charge  1 
Naval  Weather  Research  Facility 

Naval  Air  Station,  Building  R-48 
Norfolk,  Virginia  23511 

7.  Commanding  Officer  1 
U.  S.  Fleet  Weather  Central 

COMNAVMARIANAS,  Box  12 

FPO  San  Francisco,  California  96630 

8.  Commanding  Officer  1 
Fleet  Weather  Central 

Box  31 

FPO  New  York,  New  York  09540 

9.  Commanding  Officer  1 
Fleet  Numerical  Weather  Central 

Naval  Postgraduate  School 
Monterey,  California  93940 

10.  ARCRL  -  Research  Library  1 
L.  G.  Hanscom  Field 

Attn:   Nancy  Davis /Stop  29 
Bedford,  Massachusetts   01730 


66 


11.  Director,  Naval  Research  Laboratory 
Attn:  Tech.  Services  Info.  Officer 
Washington,  D.  C.   20390 

12.  American  Meteorological  Society 
45  Beacon  Street 

Boston,  Massachusetts  02128 

13.  Department  of  Meteorology 
Code  51 

Naval  Postgraduate  School 
Monterey,  California  93940 

14.  Department  of  Oceanography 
Code  58 

Naval  Postgraduate  School 
Monterey,  California  93940 

15.  Office  of  Naval  Research 
Department  of  the  Navy 
Washington,  D.  C.   20360 

16.  Commander,  Air  Weather  Service 
Military  Airlift  Command 

U.  S.  Air  Force 

Scott  Air  Force  Base,  Illinois  62226 

17.  Atmospheric  Sciences  Library 

National  Oceanographic  Atmospheric  Administration 
Silver  Spring,  Maryland  20910 

18.  Professor  Victor  Starr 
Department  of  Meteorology 
M.  I.  T. 

Cambridge,  Massachusetts  03139 

19.  Dr.  J.  Pedlosky 

Department  of  Geophysical  Sciences 
University  of  Chicago 
Chicago,  Illinois  60637 

20.  Dr.  Joanne  Simpson 
Experimental  Meteorology  Branch 

National  Oceanographic  Atmospheric  Administration 
Coral  Gables,  Florida  33124 

21.  Dr.  V.  Jurcec 

United  Nations  Development  Program 

Box  982 

Cairo,  United  Arab  Republic 


67 


22.  Dr.  A.  Huss 
Department  of  Meteorology 
Hebrew  University 
Jerusalem,  Israel 

23.  National  Center  for  Atmospheric  Research 
Box  1470 

Boulder,  Colorado  80302 

24.  Dr.  T.  N.  Krishnamurti 
Department  of  Meteorology 
Florida  State  University 
Tallahassee,  Florida  32306 

25.  Dr.  Fred  Shuman 
Director 

National  Meteorological  Center 

Environmental  Science  Services  Administration 

Suit land,  Maryland  20390 

26.  Dr.  J.  Smagorinsky 
Director 

Geophysical  Fluid  Dynamics  Laboratory 
Princeton  University 
Princeton,  New  Jersey  08540 

27.  Professor  N.  A.  Phillips 
54-1422 

M.  I.  T. 

Cambridge,  Massachusetts  02139 

28.  Dr.  E.  N.  Lorenz 
Department  of  Meteorology 
M.  I.  T. 

Cambridge,  Massachusetts  02139 

29.  Professor  J.  G.  Charney 
54-1424 

M.  I.  T. 

Cambridge,  Massachusetts  02139 

30.  Professor  K.  Ooyama 
Department  of  Meteorology 
New  York  University 
University  Heights 

New  York,  New  York  10453 

31.  Dr.  M.  G.  Wurtele 
Department  of  Meteorology 
UCLA 

Los  Angeles,  California  90024 


68 


32.  Dr.  A.  Arakawa 
Department  of  Meteorology 
UCLA 

Los  Angeles,  California  90024 

33.  Dr.  David  Houghton 
Department  of  Meteorology 
University  of  Wisconsin 
Madison,  Wisconsin  53706 

34.  Dr.  S.  K.  Kao 
Department  of  Meteorology 
University  of  Utah 

Salt  Lake  City,  Utah  84112 

35.  Dr.  J.  Holton 

Department  of  Atmospheric  Sciences 
University  of  Washington 
Seattle,  Washington  98105 

36.  Dr.  P.  Gilman 

National  Center  for  Atmospheric  Research 

Box  1470 

Boulder,  Colorado  80302 

37.  Dr.  P.  Thompson 

National  Center  for  Atmospheric  Research 

Box  1470 

Boulder,  Colorado  80302 

38.  Dr.  Peter  Stone 
Pierce  Hall 
Harvard  University 
Cambridge,  Massachusetts  02138 

39.  Dr.  John  Young 
Department  of  Meteorology 
University  of  Wisconson 
Madison,  Wisconsin  53706 

40.  Dr.  George  J.  Haltiner 

Chairman,  Department  of  Meteorology 
Naval  Postgraduate  School 
Monterey,  California  93940 

41.  Dr.  Jerry  D.  Mahlman 

Geophysical  Fluid  Dynamics  Laboratory 
Princeton  University 
Princeton,  New  Jersey  08540 


69 


42.  Dr.  Russell  Elsberry 
Department  of  Meteorology 
Naval  Postgraduate  School 
Monterey,  California  93940 

43.  Commanding  Officer 
Pacific  Missile  Range 
Attn:   Geophysics  Division 
Point  Mugu,  California  93041 

44.  Dr.  Neil  Frank 
National  Hurricane  Center 

National  Oceanographic  Atmospheric  Administration 
Coral  Gables,  Florida  33124 

45.  Lieutenant  Commander  Frank  H.  Taylor 
U.  S.  Fleet  Weather  Central 

Box  113 

FPO  San  Francisco,  California  96610 

46.  Mr.  Donald  H.  Lucas 
Fernbank  Science  Center 
156  Heaton  Park  Dr.  N.E. 
Atlanta,  Georgia  30307 

47.  Lieutenant  Commander  W.  Roger  Lambertson 
Department  of  Meteorology 

Naval  Postgraduate  School 
Monterey,  California  93940 

48.  Dr.  Peter  J.  Gierasch 
Geophysical  Fluid  Dynamics  Institute 
Florida  State  University 
Tallahassee,  Florida  32306 

49.  Dr.  Andrew  P.  Ingersoll 

Division  of  Geological  and  Planetary  Sciences 
California  Institute  of  Technology 
Pasadena,  California  91109 

50.  Dr.  Robert  L.  Haney 
Department  of  Meteorology 
Naval  Postgraduate  School 
Monterey,  California  93940 

51.  Dr.  Ron  L.  Alberty 
Department  of  Meteorology 
Naval  Postgraduate  School 
Monterey,  California  93940 


70 


52.  Dr.  W.  L.  Gates 
Department  of  Meteorology 
Naval  Postgraduate  School 
Monterey,  California  93940 

53.  Dr.  Richard  Alexander 
The  Rand  Corporation 
1700  Main  Street 

Santa  Monica,  California  90406 

54.  Commanding  Officer 
Fleet  Weather  Central 
Box  110 

FPO  San  Francisco,  California  96610 


71 


Security  Classification 


(Security   classili 


DOCUMENT  CONTROL  DATA  -R&D 

>/  title,    bodv  ol  abstract  and  indexing  annotation  must  be  entered  when   tlie 


ill  report  is   classified) 


HGINATING    ACTIVITY   ( Corporate  author) 

Naval  Postgraduate  School 
Monterey,  California  93940 


2«.   REPORT    SECURITY    C  L  A  S  SI  F  I  C  A  T  I  O  r 

Unclassified 


26.    CROUP 


REPORT     TITLE 


General  Circulation  Experiments  with  a  Two-Level 
Quasi-Geostrophic  Model  Including  Nonlinear 
Interaction  between  Two  Waves  and  the  Mean  Flow 


DESCRIPTIVE   NOTES  (Type  ol  report  andjnclusive  dates) 

Master's  Thesis;  March  1972 


S-   AUTHOR(S)  (First  name,  middle  initial,  last  name) 

Glenn  Curtis  Trumbower 


6-    REPORT    DATE 

March  1972 


7a.    TOTAL    NO.    OF    PAGES 

73 


7b.    NO.    OF    REFS 

19 


««.     CONTRACT    OR    GRANT    NO. 


b.    PROJECT    NO. 


9«.    ORIGINATOR'S    REPORT    NUMBER(S) 


10.    DISTRIBUTION    STATEMENT 


Approved  for  public  release;  distribution  unlimited. 


I  1  .    SUPPLEMEN  T  I 


12.    SPONSORING 


Naval  Postgraduate  School 
Monterey,  California  93940 


3.  ASSTR  AC  T 


Long  term  numerical   integrations  were  performed  with  a  two-level  model 
utilizing  the  quasi-geostrophic  equation  set.   Friction  was  incorporated  through 
the  vertical  derivative  of  the  eddy  stress  and  heating  was  applied  as  a  linear 
function  of  latitude.   The  long-term  interactions  between  one  or  two  waves  in  the 
zonal  direction  and  the  mean  flow  were  examined. 

The  stability  of  the  initial  mean  flow  was  investigated  with  linearized  equa- 
tions and  the  most  unstable  wave  numbers  were  determined.   The  nonlinear  equations 
were  integrated  using  the  most  unstable  wave  in  the  initial  conditions.   300  day 
forecasts  were  made  with  various  values  of  the  heating,  (3,  wall  separation,  and 
wave  number.   A  second  wave  number  was  introduced  into  the  experiments  in  order  to 
determine  the  most  likely  wave  number  for  the  fully  evolved  mean  flow.   Constant 
amplitude,  baroclinic,  propagating  disturbances  were  obtained  in  every  case  where 
the  wall  separation  was  4,000  km.   The  experiments  with  the  wall  separation  of 
8,000  km  produced  non-steady  waves  with  fluctuating  baroclinic  and  barotropic 
interactions. 


DD,F°r.,1473 

S/N    01 01 -807-681 1 


(PAGE     1  ) 


72 


Security  Classification 


Security  Classification 


EY  WORDS 


baro clinic 

baro tropic 

constant  amplitude 

eddy  stress 

energy  transformations 

general  circulation 

quasi-geos trophic 

two-level  model 

vorticity 


DD,F°oR.\,1473  "»«) 

S/N     U1 01 -807-682 1 


73 


Security  Classification 


T8223 
c.l 


Trumbower 

Genera!  circulation 
experiments  with  a  two- 
level  Quasi -Geos trophic 
model  including  non- 
linear interaction  be- 
tween two  waves  and  the 
mean  flow. 


Thesis  .1.34943 

T3223  Trumbower 

c.l        General  circulation 
experiments  with  a  two 
level  Quasi -Geos trophic 
model  including  non- 
linear interaction  be- 
tween two  waves  and  the 
mean  flow. 


thesT8223 

General  circulation  experiments  with  a  t 


3  2768  001  88847  2 

DUDLEY  KNOX  LIBRARY 


